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Abstract 

Feedback effects due to massive stars and supernovae 
in the first objects are shown to strongly regulate both 
galaxy formation/evolution and the reionization process. 
Here we review the most important ones in some detail. 
We discuss how Type II supernovae can be used as tracers 
of the first objects and detected with NGST, for which we 
predict supernova number counts including the effects of 
gravitational lensing. Preliminary results on the formation 
of dust in the ejecta of supernovae of primordial compo- 
sition are also presented. We finally turn to the consid- 
eration of the process of inhomogeneous reionization due 
to primordial stellar sources by means of high resolution 
numerical simulations, allowing for a self-consistent treat- 
ment of the above feedback processes. These simulations 
allow us to draw conclusions on the evolution and epoch 
of reionization and about the fate of reionizing objects. 
We conclude that a large fraction (« 99%) of collapsed 
objects must be dark at redshift around eight. 



1. The First Objects 

As the temperature of the cosmic bath decreases, atoms 
start to recombine and therefore decouple from CMB ra- 
diation at redshift « 1100. The baryonic Jeans mass after 
this event is given by (assuming f2 = 1) 

3/2 / rp \ 3/2 

j 1],Mq, (1) 



Mj ~ 6 X 10 
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30 J V500K, 
where T is the gas temperature and the baryon den- 
sity parameter. Masses larger than Mj are gravitationally 
unstable and should, in principle, collapse. However, in 
order for the actual collapse to occur a more severe con- 
dition must be satisfied, i.e. that the cooling time of the 
gas is shorter than the Hubble time at that epoch. In fact, 
radiative losses provide the only way for the gas to lose 
pressure and to settle down in the potential well of the 
dark matter halo. Since the virial temperature correspond- 
ing to the masses of the first objects (PopIII) is typically 
8000 K, cooling by hydrogen Lya excitation is strongly 
quenched, and the only viable coolant in a primordial H- 
He plasma is molecular hydrogen. H2 is produced during 
the recombination phase, but its relic abundance is very 
small. Primordial H2 forms with a fractional abundance of 



« 10""^ at redshifts > 400 via the H^ formation channel. 
At redshifts ^ 110, when the Cosmic Microwave Back- 
ground radiation (CMB) intensity becomes weak enough 
to allow for significant formation of H^ ions, a primordial 
fraction of /^^ ~ 2x 10"^ (GaUi & Palla 1998) is produced 
for model universes that satisfy the standard primordial 
nucleosynthesis constrain ^i,h? = 0.0125 (Copi, Schramm 
& Turner 1995), where fib is the baryon density parame- 
ter and i/o = lOO/i km s^^ Mpc^^ is the Hubble constant. 
This primordial fraction is usually lower than the one re- 
quired for the formation of Pop III objects, but during the 
collapse phase the molecular hydrogen content can reach 
high enough values to trigger star formation. On the other 
hand, objects with virial temperatures (or masses) above 
that required for the hydrogen Lya line cooling to be ef- 
ficient, do not rely on H2 cooling to ignite internal star 
formation. Thus, the the fate of a virialized lump depends 
crucially on its ability to rapidly increase its H2 content 
during the collapse phase. Tegmark et al. (1997) have ad- 
dressed this question in great detail by calculating the evo- 
lution of the H2 abundance for different halo masses and 
initial conditions for a standard CDM cosmology. They 
conclude that if the prevailing conditions are such that a 
molecular hydrogen fraction of order of //j^ « 5 x 10^^ 
is produced during the collapse, then the lump will cool, 
fragment and eventually form stars. This criterion is met 
only by larger halos implying that for each virialization 
redshift there will exist some critical mass. Merit , such 
that protogalaxies with total mass M > Merit will be 
able to form stars and those with M < Merit will fail. 
As an example, a 3a fluctuation of a Cold Dark Matter 
primordial spectrum, has Merit ~ IO^Mq and collapses at 
z « 30. 

The first objects affect the subsequent galaxy forma- 
tion, reionization and metal enrichment of the universe 
through their feedback effects. Generally speaking, two 
types of feedback, radiative and stellar, can be at work. 
They are related to the ionizing/ H2 -photodissociating 
radiation and mechanical energy input from massive stars 
and supernovae, respectively. We discuss them in detail in 
the following. 

2. Radiative Feedback 

As stars are formed inside the first collapsed objects, their 
ionizing radiation produces HII regions in the surround- 
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ing IGM, thus triggering the reionization of the universe. 
In addition, their UV photons start to photodissociate 
H2 molecules in the neighbor objects. We have seen that 
H2 is a crucial species for the cooling to occur, and the 
lack of it might prevent the collapse of small protogalaxies. 
This process is known as radiative feedback. Particularly 
relevant is the soft UV background in the LW bands, as by 
dissociating the H2, it could influence the star formation 
history of other small objects preventing their cooling. 

Ciardi, Ferrara & Abel (2000) have shown that the 
UV flux from these objects results in a soft (Lyman - 
Werner band) UV background (SUVB), Jlw, whose in- 
tensity (and hence radiative feedback efflciency) depends 
on redshift. At high redshift the radiative feedback can be 
induced also by the direct dissociating flux from a nearby 
object. In practice, two different situations can occur: i) 
the collapsing object is outside the dissociated spheres 
produced by preexistent objects: then its formation could 
be affected only by the SUVB {JLW,b), as by construction 
the direct flux {JLW,d) can only dissociate molecular hy- 
drogen on time scales shorter than the Hubble time inside 
this region; ii) the collapsing object is located inside the 
dissociation sphere of a previously collapsed object: the 
actual dissociating flux in this case is essentially given by 
JLW,max = {J Lw,b + J Lw,d) ■ ^ t^us assumed that, given 
a forming Pop III, if the incident dissociating flux {JLW,b 
in the former case, JLW.max in the latter) is higher than 
the minimum flux required for negative feedback (J^), the 
collapse of the object is halted. This implies the existence 
of a population of "dark objects" which were not able to 
produce stars and, hence, light. 

To assess the minimum flux required at each redshift 
to drive the radiative feedback Ciardi et al. (2000, CFGJ) 
have performed non-equilibrium multifrequency radiative 
transfer calculations for a stellar spectrum (assuming a 
metallicity Z = 10^''^) impinging onto a homogeneous 
gas layer, and studied the evolution of the following nine 
species: H, H~, 11+ , He, He+, He+"'", H2, Hj and free 
electrons for a free fall time. We can then define the mini- 
mum total mass required for an object to self-shield from 
an external flux of intensity J^fl at the Lyman limit, as 
Msh = (4/3) i:{ph)Rlt^, where {pn) = {pb)^b^ is the mean 
dark matter density of the halo in which the gas collapses; 
Rsh is the shielding radius beyond which molecular hydro- 
gen is not photodissociated and allowing the collapse to 
take place on a free-fall time scale. 

Values of Msh for different values of Js^ have been 
obtained at various redshifts (see Fig. |l|) and will be then 
used in the reionization calculations discussed below. Pro- 
togalaxies with masses above Mh for which cooling is 
predominantly contributed by Lya line cooling will not 
be affected by the negative feedback studied here. The 
collapse of very small objects with mass < Merit is on 
the other hand made impossible by the cooling time be- 
ing longer than the Hubble time. Thus the only mass 
range in which negative feedback is important lies ap- 



proximately in 10^ — 1O*M0, depending on redshift. In 
order for the negative feedback to be effective, fluxes of 
the order of 10"^'' — lO^^'^erg s^^ cm^^ Hz^^ sr"^ are 
required. These fluxes are typically produced by a Pop III 
with baryonic mass IO^M^^sMq at distances closer than 

1 /2 

~ 21 — 7 X kpc for the two above flux values, re- 

spectively, while the SUVB can reach an intensity in the 
above range only after z « 15. This suggests that at high z 
negative feedback is driven primarily by the direct irradia- 
tion from neighbor objects in regions of intense clustering, 
while only for z < 15 the SUVB becomes dominant. 




10-2S 10-24 IQ-Z2 

Jso [^^S s ' cm-2 Hz"' sr^'] 



Figure 1. Minimum total mass for self- shielding from an 
external incident flux with intensity Js,o o-t the Lyman 
limit. The curves are for different redshift: from the top 
to the bottom z =8.3, 11. 4, 15.4, 18.0, 22.1, 29.6. Circles 
show the value of Mu (open) and Merit (filled). Radiative 
feedback works only in the solid portions of the curves. 



3. Stellar Feedback 

Once the first stars have formed in the host protogalaxy, 
they can deeply influence the subsequent star formation 
process through the effects of mass and energy deposi- 
tion due to winds and supernova explosions. While low 
mass objects may experience a blowaway, expelling their 
entire gas content into the IGM and quenching star for- 
mation, larger objects may instead be able to at least par- 
tially retain their baryons. However, even in this case the 
blowout induces a decrease of the star formation rate due 
to the global heating and loss of the galactic ISM. These 
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two regimes are separated by a critical mass, Mi,y (CFGJ, 
Ferrara & Tolstoy 2000, MacLow & Ferrara 1999). For 
the relatively small objects present during the reionization 
epoch 30 ^ ^ 10 the importance of these stellar feed- 
hacks can hardly be overlooked. To understand the role 
of stellar feedback, let us consider a collapsed object with 
mass lower than M^y. Then the star formation is suddenly 
halted as the entire gas content is removed. In this case, 
the ionizing photon production will last only for a time 
interval of order toB ~ 10^ yr^ the mean lifetime of the 
massive stars produced initially. After this phase, the ion- 
ized gas around the source will start to recombine at a 
fast rate as a result of the highly efficient high z Compton 
cooling, rapidly decreasing the temperature inside the ion- 
ized region. Due to the short lifetime and recombination 
time scales, these object will only produce transient HII 
regions which will rapidly disappear. In fact, the recom- 
bination time scale, when the Compton cooling is taken 
into account, is of order tree « (1 — 50) x 10^ yr at redshift 
2:«30-10, respectively, and therefore much shorter than 
the corresponding Hubble time. The dissociating photon 
production will nevertheless continue for a longer time, 
due to the important contribution of long-lived interme- 
diate mass stars formed in the same initial star formation 
burst. This, combined with the fact that there is no effi- 
cient mechanism available to re-form the destroyed H2 in 
the IGM analogous to H recombination, implies that the 
dissociation is not impeded by blowaway and the contribu- 
tion from these small objects should be yet accounted for. 
However, this is not necessary. In fact, after blowaway, H2 
is efficiently formed in the shocked IGM gas, cooling un- 
der non equilibrium conditions (Ferrara 1998). The final 
radius of the cooled shell behind which molecular hydro- 
gen is formed by this process is: 



-19/10 



kpc. 



(2) 



We notice that the formation process produces an amount 
of H2 roughly similar to the one destroyed by radiation 
(Ferrara 1998). A lower limit to the amount of molecular 
hydrogen produced in an explosion is readily found to be 
equal to M+^{z) ~ (4^/3)pfci?3(2/) or 

M+Jz) ^ lQ2fl'in,,,{l + z)-^ /^"il^/Ve/i'/'^Me. (3) 

The UV/ionizing radiation from Pop HI massive stars pre- 
vious to blow-away will produce both a HII region and 
a region of photodissociated intergalactic H2 (radius Rd) 
in which the object is embedded. The radius Rd can be 
defined by requiring that the photodissociation timescale 
td is shorter than t^- This condition yields the defini- 
tion Rd = S]I^{1 + z)-3/4/j-i/2^ ^i^ere Slw = PSi{{)) is 
the UV photon flux in the H2 Lyman- Werner (LW) bands 
(11.2-13.6 eV), assumed here to be proportional to the 
flux of Lyc photons, S'i(O), just before the massive stars 
explode. The value of the constant (3 depends somewhat 
on the IMF and on the evolutionary stage of the stellar 



cluster, but its value should be reasonably close to unity. 
We estimate S'i(O) to be 

S'i(O) ~ 10*'^ fiivpp,48fesc,20^b,5fb,sM5 S~"^, (4) 

where fuvpp = 48000/„t,pp_48 is the UV photon produc- 
tion per collapsed proton efficiency (Tegmark et al. 1994) 
and fesc = 0.2/esc,20 is the escape fraction of such pho- 
tons. It is worth noticing that the latter quantity is rather 
uncertain. Recently, Dove, ShuU & Ferrara (2000) have 
estimated the fraction of ionizing photons emitted by OB 
associations that escapes the H I disk of our Galaxy into 
the halo and intergalactic medium (IGM) by solving the 
time-dependent radiation transfer problem of stellar ra- 
diation through evolving superbubbles within a smoothly 
varying H I distribution. They find that the shells of the 
expanding superbubbles quickly trap or attenuate the ion- 
izing flux, so that most of the escaping radiation escapes 
shortly after the formation of the superbubble. Superbub- 
bles of large associations can blow out of the H I disk 
and form dynamic chimneys, which allow the ionizing ra- 
diation to escape the H I disk directly. However, blowout 
occurs when the ionizing photon luminosity has dropped 
well below the association's maximum luminosity. For a 
coeval star-formation history, the total fraction of Lyman 
Continuum photons that escape both sides of the disk in 
the solar vicinity is fesc ~ 0.15 ±0.05; for a Gaussian star 
formation history, fesc ~ 0.06 ± 0.03. 

With these assumptions, it follows that 

Rd ^ iPSi,y^'il + z)^^'^h-^'^ kpc. 



^30 " "^i'^i (5) 

where 6*47 — S'i(0)/10''^s^^. Thus, the ratio between the 
H2 mass produced and destroyed by Pop Ills is 



/ 



I IGM 



Rs 



n 3 



Rd 



(6) 



The previous relation is graphically displayed in Fig. |^ 
along with the values of Rs and Rd- From that plot we see 
that objects of total mass M5 = 10 produce more H2 than 
they destroy for z ^ 25; larger objects (M5 = 100) provide 
a similar positive feedback only for 2: ^ 15, since they are 
characterized by a higher Rd/Rs ratio. However, since in a 
hierarchical model larger masses form later, even for these 
objects the overall effect should be a net H2 production. 

As inside Rs H2 is re-formed, the net effect on the 
destruction of intergalactic H2 of blown away objects is 
negligible (if not positive) and we can safely neglect them 
in the subsequent calculations. 

In conclusion, low mass objects with mass < Mfjy pro- 
duce ionization regions which last only for a recombina- 
tion time; larger objects can survive but their star forma- 
tion ability is impaired by ISM loss/heating. The detailed 
derivation of Mby can be found in CFGJ. 

4. Evolutionary Tracks due to Feedbacks 

Fig. ^ illustrates all possible evolutionary tracks and final 
fates of primordial objects, together with the mass scales 
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Figure 2. Ratio between the H2 mass formed and destroyed 
by a Pop III object as a function of the SNII explosion red- 
shift (solid lines); values larger than one for the ratio de- 
fine the epochs where Pop Ills have a positive feedback on 
galaxy formation. Also shown are the shell (proper) radius 
at cooling, Rg (dotted), and the photodissociation (proper) 
radius, Rd (dashed). The upper set of curves refers to ob- 
jects of mass M — 10^ Mq, whereas the bottom one corre- 
sponds to larger objects, M = lO^M©. 



determined by the various physical processes and feed- 
backs. We recall that there are four critical mass scales in 
the problem: (i) Merit, the minimum mass for an object 
to be able to cool in a Hubble time; (ii) Mh, the critical 
mass for which hydrogen Lya line cooling is dominant; 
(iii) Msh, the characteristic mass above which the object is 
self-shielded, and (iv) M^y the characteristic mass for stel- 
lar feedback, below which blowaway can not be avoided. 
Starting from a virialized dark matter halo, condition (i) 
produces the first branching, and objects failing to sat- 
isfy it will not collapse and form only a negligible amount 
of stars. In the following, we will refer to these objects 
as dark objects. Protogalaxies with masses in the range 
Merit < M < Mh are then subject to the effect of ra- 
diative feedback, which could either impede the collapse 
of those of them with mass M < Msh, thus contributing 
to the class of dark objects, or allow the collapse of the 
remaining ones (M > Mgh) to join those with M > Mh 
in the class of luminous objects. This is the class of ob- 
jects that convert a considerable fraction of their baryons 



Figure 3. Possible evolutionary tracks of objects as deter- 
mined by the processes and feedbacks discussed in this pa- 
per 

in stars. Stellar feedback causes the final bifurcation by 
inducing a blowaway of the baryons contained in lumi- 
nous objects with mass M < Mby-, this separates the class 
in two subclasses, namely "normal" galaxies (although of 
masses comparable to present day dwarfs) that we dub 
gaseous galaxies and tiny stellar aggregates with negligi- 
ble traces (if any) of gas to which consequently we will 
refer to as naked stellar clusters. 

5. Detecting PopIII Objects 

Before we proceed to show how the reionization of the uni- 
verse is regulated by stellar feedbacks, we briefly discuss 
possible strategies to observe PopIII objects. Given their 
small mass, these objects are likely to be faint: for a rea- 
sonable mass-to-light ratio for young galaxies « 0.1, their 
bolometric luminosity is w 2 x 10"*" erg s~^. A Type II SN 
is typically one hundred times brighter. Thus, for peri- 
ods even longer than a year (taking into account the time 
stretching c>c (1 + z) of the light curve) the SNII outshines 
its host galaxy. Since Type II SNe originate from massive 
stars, a necessary requirement is that the (unknown) IMF 
of these objects is flat enough to extend into this regime. 
Clearly, this condition is satisfied in the local universe, but 
it is not necessarily so in the conditions prevailing when 
the universe was young. Nevertheless, it is intriguing to 
speculate about the observational perspectives to detect 
very high z SNe, a possibility to which our hopes to in- 



5 



vestigate directly the primeval star/galaxy formation are 
closely tied. 

The advances in technology are making available a 
new generation of instruments, some already at work and 
some in an advanced design phase, which will dramatically 
increase our observational capabilities. As representative 
of such class, we will focus on a particular instrument, 
namely the Next Generation Space Telescope (NGST) . In 
the following we will try to quantify the expectations for 
the detection of high z SNe. In addition to what we can 
learn on the formation of the first objects. Pop III SNII 
can provide crucial information on cosmological models 
due to their different predictions concerning the gravita- 
tional lensing magnification patterns by the intervening 
matter in the universe. The gravitational amplification 
might affect the SNII number counts and should be there- 
fore properly taken into account (Marri & Ferrara 1998). 
As an example, taken from the study of Marri, Ferrara & 
Pozzetti (2000), we show the magnification map for the 
SCDM model for a source located at redshift 10 in Fig. 
^. Fig. H shows the differential SNII counts [0.5 mag/ yr/ 




Figure J^. Magnification map for a SCDM model 4\ 
corresponding to a NGST field) for a source located at 
z = 10. The magnification range is 0.7-50. 



0.44 deg^] as a function of AB magnitude. The four panels 
contain the curves for the SCDM and LCDM models and 
for J, K, L, and M bands, both including or neglecting the 
effects of GL. For comparison, we plotted the NGST mag- 
nitude limit AB = 31.4 (vertical line). This is calculated 
by assuming a constant limiting flux Tngst = 10 nJy in 
the wavelength range 1 — 5/xm {i.e. J—M bands). This can 



be achieved, for a 8-m (10-m) mirror size and a S/N=5, 
in about 2.6 x 10^ s (1.1 x 10^ s)[]. Thus, NGST should 
be able to reach the peak of expected SNII count dis- 
tribution, which is located at AS « 30 - 31 for SCDM 
and AB k, 31 — 32 for LCDM (depending on the wave- 
length band). The differences among the various bands 
are not particularly pronounced, although J and K bands 
present a larger number of luminous {AB ^ 27) sources, 
and therefore they might be more suitable for the exper- 
iment. Furthermore, we point out that in the L and M 
bands NGST, with the current magnitude limit, will not 
be able to reach the peak of expected SNII count distri- 
bution in the LCDM model. 




Figure 5. Differential number counts for the two cosmolog- 
ical models considered (thick solid line: SCDM, thin solid 
line: LCDM) as a function of apparent AB magnitude in 
J ,K ,L,M hands. Also shown is the NCST limiting mag- 
nitude. Dashed curves neglect lensing magnification. 



6. More feedback from dust formation 

In addition to the above mentioned feedback effects, SNII 
can also initiate molecular hydrogen formation on dust 
grain surfaces rather than in the gas phase, a more stan- 
dard process in present day galactic environments. At high 
redshift Type II SNe are the only possible sources of dust, 
due to the short age of the universe and the long evolu- 
tionary timescales characterizing more conventional dust 
sources, as for example evolved stars. A detailed nucle- 
ation study of dust formation in the ejecta of Type II SNe 



^ This result has been obta ined using the NGST Ex posure 
Time Calculator, available at jittp:/ /augusta.stsci.edu 
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whose initial chemical composition is primordial has been 
very recently carried on (Todini & Ferrara, in prepara- 
tion). This case is suitable for the study of the very first 
SNII, as those originating in PopIII objects. From this 
study we have been able to derive several important prop- 
erties characterizing the first solid particles produced in 
the universe, as for example the total dust mass produced, 
the dust composition and the grain size distribution. Pre- 
liminary results show that in the range of SNII masses 
considered, Msat = 12 - 40Mq, typically 0.1 Mg of dust 
is produced. SNII towards the low end of the mass distri- 
bution mostly produce silicate grains, whereas the most 
massive ones predominantly produce graphite grains; this 
fact is easy to understand in terms of the C/0 ratio of 
their ejecta. The size distribution is approximated by a 
power-law for about two decades in radius (from « lA to 
lOOA) and shows a smooth cutoff beyond that grain ra- 
dius. Thanks to these results we can now pose the fol- 
lowing question: what is the minimum amount of dust 
required in order for the molecular hydrogen formation 
on grains to become competitive with the gas phase one 
? An order-of-magnitude answer can be obtained by com- 
paring the two formation rates. At the low densities rel- 
evant here, II2 is formed in the gas phase mainly via the 
channel H + e~ ^ H~ + hi/, at rate fcg (the rate coefh- 
cient fcg is given in Abel et al. 1997); formation via the 
H2 ^ channel, when included, is found to be negligible in 
our case. Therefore the formation rate in the gas phase 
is 7?. ~ ksrifj-nH- The formation rate on grain surfaces 
is instead given by TZ"^ ~ 0.5{'~fCs(j)ndnH , where 7 is the 
sticking coefficient, Cg is the sound speed in the gas, and 
a is the grain cross section. The equality between the two 
rates can be cast into the following form: V = O.l^/Txe, 
where V is the dust-to-gas ratio normalized to its Galac- 
tic value, T is the gas temperature and Xg the gas ioniza- 
tion fraction. For typical parameters of the PopIII objects, 
H2 production on dust grains becomes dominant once V is 
larger than 5% of the local value. With the dust yields per 
SNII calculated above we then conclude that only about 
15 SNII are required to enrich in dust to this level a pri- 
mordial object. Clearly, early dust formation might play 
a role in the formation of the first generation of objects. 

7. Reionization History Driven by Feedbacks 

In this final Section we summarize some of the results ob- 
tained by CFGJ on the evolution of inhomogeneous reion- 
ization, regulated by the feedback mechanisms introduced 
above and regulated by the molecular hydrogen network. 
For a broader presentation the reader is advised to re- 
fer to CFGJ. We have determined the spatial distribu- 
tion of the ionizing sources from high resolution numeri- 
cal N-body simulations within a periodic cube of comoving 
length L = 2.55h~^ Mpc for a critical density cold dark 
matter model (r2o = li h=0.5 with (78=0.6 at z=0). This 
allows us to describe the topology of the ionized and dis- 



sociated regions at various cosmic epochs and derive the 
evolution of H, He, and H2 filling factors, soft UV back- 
ground, cosmic star formation rate and the final fate of 
ionizing objects. There are three free parameters in the 
computation: /(,, the fraction of virialized baryons that is 
able to cool and become available to form stars, /* the 
star formation efficiency, and fesc the photon escape frac- 
tion from the proto-galaxy. The variation of these free pa- 
rameters defines four different runs A-D, whose parameter 
combination is shown in Tab. 1. 



RUN 


fb 


,A 


fesc 


A 


0.08 


0.15 


0.2 


B 


0.08 


0.05 


0.2 


C 


1.00 


0.15 


0.2 


D 


0.08 


0.15 


0.1 



Table 1. Parameters of the calculation: fraction of virial- 
ized baryons that are able to cool and be available for star 
formation, fb; star formation efficiency, ft,; photon escape 
fraction, fesc- 



7.1. Filling factors 

The filling factors of the dissociated H2 and ionized H 
and He are defined as the box volume fraction occupied 
by those species. The results are shown in Fig. ^ for the 
different runs. The intergalactic relic molecular hydrogen 
is found to be completely dissociated at very high redshift 
[z Ri 25) independently of the parameters of the simula- 
tion. This descends from the fact that dissociation spheres 
are relatively large and overlap at early times. Ionization 
spheres are instead always smaller than dissociation ones 
and complete reionization occurs considerably later. Ex- 
cept for run C, when reionization occurs by z~15, pri- 
mordial galaxies are able to reionize the IGM at a redshift 
z^lO. The filling factor is approximately linearly propor- 
tional to the number of sources but has a stronger (cu- 
bic) dependence on the ionization sphere radius. Thus, 
although the number of (relatively small) luminous ob- 
jects present increases with redshift, the filling factor is 
only boosted by the appearance of larger, and therefore 
more luminous, objects which can ionize more efficiently. 
The subsequent flattening at lower redshifts, is obviously 
due to the fact that when the volume fraction occupied 
by the ionized gas becomes close to unity most photons 
are preferentially used to sustain the reached ionization 
level rather than to create new ionized regions. The reion- 
ization is basically driven by objects collapsed through H 
line cooling, while small mass PopIII stars play only a 
minor role and even in the absence of a radiative neg- 
ative feedback they would not be able to reionize the 
IGM. Note that only about 2% of the stars observed at 
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Figure 6. Evolution of the dissociated molecular hydrogen 
(upper set of lines), ionized atomic hydrogen (middle set) 
and doubly ionized helium (bottom line) filling factor as 
a function of redshift for different runs: A (solid line), B 
(dotted), C (dashed- dotted) and D (dashed), (b) same as 
(a) for runs: A (solid line), Al (dotted) and A2 (dashed- 
dotted). 

z = is required to reionize the universe completely. This 
corresponds to an average IGM metallicity (assuming a 
Salpeter IMF) at redshift z«10 equal to (Z) w 3x IQ-'^Zq. 

7.2. What is the fate of (re)ionizing sources ? 

The last point we make here concerns the final fate of 
the objects responsible for the reionization of the uni- 
verse in terms of the evolutionary tracks introduced above; 
this discussion is summarized in Fig. |^. The straight lines 
represent, from the top to the bottom, the number of 
dark matter halos, dark objects, naked stellar clusters and 
gaseous galaxies, respectively. The dotted curve represents 
the number of luminous objects with large enough mass 
(M > Mh) to make the H line cooling efficient and be- 
come insensitive to the negative feedback. We remind that 
the naked stellar clusters are the luminous objects with 
M < Mby, while the gaseous galaxies are the ones with 
M > Mby; thus, the number of luminous objects present 
at a certain redshift is given by the sum of naked stel- 
lar clusters and gaseous galaxies. We first notice that the 
majority of the luminous objects that are able to form at 
high redshift will experience blowaway, becoming naked 
stellar clusters, while only a minor fraction, and only at 
z ^ 15, when larger objects start to form, will survive and 
become gaseous galaxies. An always increasing number of 
luminous objects is forming with decreasing redshift, until 



Figure 7. Number evolution of different objects in the sim- 
ulation box for run A. 

z«15, where a flattening is seen. This is due to the fact 
that the dark matter halo mass function is still dominated 
by small mass objects, but a large fraction of them can- 
not form due to the following combined effects: i) towards 
lower redshift the critical mass for the collapse (Merit) in- 
creases and fewer objects satisfy the condition M > Merit', 
ii) the radiative feedback due to either the direct dissoci- 
ating flux or the SUVB increases at low redshift as the 
SUVB intensity reaches values significant for the negative 
feedback effect. When the number of luminous objects be- 
comes dominated by objects with M > Mh, by zwlO the 
population of luminous objects grows again, basically be- 
cause their formation is now unaffected by negative feed- 
back. A steadily increasing number of objects is prevented 
from forming stars and remains dark; this population is 
about «99% of the total population of dark matter halos 
at z«8. This is also due to the combined effect of points i) 
and ii) mentioned above. This population of halos which 
have failed to produce stars could be identified with the 
low mass tail distribution of the dark galaxies that reveal 
their presence through gravitational lensing of quasars. It 
has been argued in the recent literature that this popu- 
lation of dark galaxies outnumbers normal galaxies by a 
substantial amount, and Fig. M supports this view. 



8. Summary 

We have discussed several feedback effects, mostly related 
to the presence of massive stars in early formed objects, 
and analyzed their impact on galaxy formation and reion- 
ization of the universe. A brief summary of the findings of 
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this paper (for a more detailed presentation see Ciardi et 
al. 2000) follows below: 

• The first objects affect subsequent galaxy formation 
both via radiative and stellar feedbacks due to ionizing/ 
H2 -photodissociating radiation and mechanical energy 
input due to massive stars and supernovae, respectively 

• Molecular hydrogen is photodissociated but reformed 
behind SN shocks during blowaway/blowout events; dust 
can be formed in SNII ejecta, allowing for H2 formation 
on grain surfaces rather than in the gas phase. 

• PopIII objects could be detected via their SNII with 
NGST; Up to z = 15, the (SCDM, LCDM) models predict 
a total number of (857, 3656) SNII/yr in 100 surveyed 
4' X 4' fields of the Next Generation Space Telescope. 

• Galaxies are able to reionize the neutral atomic hy- 
drogen by a redshift z«10, while molecular hydrogen is 
completely dissociated at very high redshift (z«25). 

• IGM reionization is basically driven by objects col- 
lapsed through H line cooling (M > Mh), while small 
mass objects {M < Mh) play only a minor role and even 
in the absence of a radiative negative feedback they would 
not be able the reionize the IGM. 

• A consistent fraction of halos is prevented from form- 
ing stars by either the condition M < Alcrit or the effect of 
radiative feedback; this population of dark objects reaches 
«99 % of the dark matter halo population at z«8. 
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